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Abstract. A sample of 54 selected Active Galactic Nuclei (AGN) has been observed with the HEGRA stereoscopic system 
of Cherenkov Telescopes between 1996 and 2002 in the TeV energy regime. The observations were motivated by the positive 
results obtained for Mkn 421 and Mkn 501. The distances of the selected objects vary over a large range of redshifts between z 
= 0.004 and z = 0.7. Among the observed AGN are the meanwhile established TeV-emitting BL Lac type objects H 1426+428 
and lES 1959+650. Furthermore the BL Lac object lES 2344+514 and the radio galaxy M 87 show evidence for a signal on 
a 4 a" level. The observation of lES 2344+514 together with the Whipple results firmly establishes this AGN as a TeV source. 
Several objects (PKS 2155-304, BL Lacertae, 3C066A) that have been claimed as TeV y-ray emitters by other groups are 
included in this data sample but could not be confirmed using data analysed here. The upper limits of several AGN included in 
this analysis are compared with predictions in the frame-work of SSC models. 

Key words. Gamma rays: observations - Galaxies: active - BL Lacertae objects: individual: lES 2344+514 - BL Lacertae 
objects: individual: lES 1959+650 - BL Lacertae objects: individual: H 1426+428 ~ radio galaxies: individual: M 87 

1. Introduction detection of TeV y-rays from the radio galaxy M87 with the 

HEGR A Cherenkov telescopes was reported ( Aharonian et alJ 
In the commonly adopted view the 'central engine' of Active |2003ij) 

Galactic Nuclei (AGN) consists of a super massive black 

hole with up to 10^ Mq surrounded by an accretion disk. 

Two relativistic plasma outflows (jets) perpendicular to the , , ■ ■ . . . r,. , , • 

,. , , , , . . _ . the IC mechanism is assumed to produce the TeV emis- 
accretion disk can be observed in some AGN ( Rees 1984; , — ' ' ' 



Different models for the production of TeV y-rays from 
BL Lac objects have been proposed. In leptonic models 



r: p „ . 1 L„„ J „ „ . ;r7r sion (e.g. 'Sikora^ J2001I) ). whereas in hadronic models the y- 

lUrrv & Padovanil 119951) . So far, y-rays in the TeV energy i-, -- 
rays are produc ed via the interactions of relativistic protons 



regime from AGN have essentially been detected from obiects ., , |_ , , . „ , ,. , — ■ — 1 1-„„^, ,. , 

°, , , ■ , • • , with matter (e.g.lPohl & Schlickeisei J2000I)). ambien t photons 

of the BL Lac type, i.e. AGN having their let pointing close rr. ; — • i, r.r.AL u \ i ■ l-.r.r.nl. i i 

.c .r^, I,, r,^,. (iMannheim 1993') or magnetic field i AharoniaDl200(J) . or both 

to the observer s line of sight. Furthermore, all known TeV ^ 

blazars are X-ray selected BL Lac objects. Recently, the first 



Mucke & Prot heroe20QU - 

The observed TeV emission shows high flux variability 
Send offprint requests to: M. Tluczykont, e-mail: on timescales stretching from months to less than an hour. 
Martin.Tluczykont@desy.de Detailed studies of variability of BL Lac type objects can 
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contribute to the understandin g of their in trinsic accelera- 
tion rn echanisms (Krawczvnski et al] l200lt [\haronian et al. 

The positive results obtained from the observa- 
tions of t he prominent extra g alactic sources of TeV y-rays 
Mkn4 21 jPunch et alJ Il992t IPetrv et all Il996t IPiron et alJ 
200 ih and Mkn501 llOuinn et alJll996t iBradburv et al.lll997t 
Dian nati-At ai et all 1 19991) as well as their relevance for 
the questio n of t he extragalactic background photon field 
jAharonianl l200ll) . have motivated further observations of 
Active Galactic Nuclei with the HEGRA Cherenkov tele- 
scopes. In this paper we present the results of dedicated obser- 
vations of 54 AGN in the years 1996 to 2002 with the stereo- 
scopic system of Cherenkov telescopes. After a brief introduc- 
tion to the HEGRA Cherenkov telescopes the analyzed data set 
will be presented followed by a description of the analysis used 
in this paper and a presentation of the results. The paper closes 
with a discussion and a summary. 

2. The HEGRA Cherenkov telescope system 

The HEGRA stereoscopic Cherenkov telescope system (1996- 
2002) co nsisted of 5 ima ging air Cherenkov telescopes 
(lACTs) dPaum et alJ Il997i) used in the stereoscopic obser- 
vation mode on the Canary island of La Palma (28.75° N, 
17.90° W) at an altitude of 2 200 m a.s.l. Additionally, one tele- 
scope (not u sed for the present an alysis) was operated in stand- 
alone mode (^ Mirzovan et al.^ 1994). Each system telescope was 
equipped with an 8.5 m^ tesselated mirror dish of 30 single mir- 
rors with a diameter of 60 cm each, and a camera consisting of 
271 photomultiplier tubes (pixels). The HEGRA lACT system 
was operating at an energy threshold of 0.5 TeV for photons 
of vertical incidence, with energy and angular resolution of 
AE/E = 10 - 20 % and 0.1° respectively on an event-by-event 
basis. The field of view of each system telescope had a diameter 
of 4.3°. The introduction of the stereoscopic observation tech- 
nique results in an improvement of the sensitivity of Cherenkov 
telescopes and especially allows for an effective y-hadron sep- 
aration (see below). The performance of the HEGRA system o f 
Cherenkov telescopes can be found in |PuhlhoferetaiTll2003h . 

3. Data Set 

Observations of 54 objects of the AGN class were carried out 
from 1996 to 2002, resulting in a total pre-selection expo- 
sure time of approx. 1150 hours (not including Mkn421 and 
Mkn501) corresponding to more than one year of continuous 
observations in moonless nights with the HEGRA lACT sys- 
tem. The total observation time accumulated for Mkn421 and 
Mkn501 amounts to more than 1500 hours. The results of the 
HEGRA observations f rom these two obiects were presented 
in different publications (Bradbury et al. '1997";' Aharonian et al. 
1999a b c: Sambruna et al. 2000: Aharonian et al.. .2001 a.b- 
2002ah and are not included in the present work. 

Mkn421 and Mkn501 excluded, the data set contains 
37 objects identified and confirmed as BL Lac type ob- 
jects, 10 radio galaxies, 4 Seyfert ga laxies, 1 qua s ar an d 
2 galaxies following the catalogues of Stickel et al. (1994); 
IPadovani & Giommi (.1995.) : .Veron-Cettv & Veron. (.2001.) . 



The object lES 0806-t-524 is one of the BL Lac objects which 
was proposed by Tinyakov & Tkachev (2001) to coincide with 
an AGASA UHECR triplet, thus being a candidate for the 
acceleration of ultra high energy cosmic rays. The distances 
of the objects vary over a large range of redshifts between 

0.004 (M87) and z*0.7 (PKS 0219-164). However, the 
expected absorption due to pair production of TeV y-rays 
with the extragalactic background light (EBL) ( lNikishovll962l) 
(TtcV Tebl — > e^e") increases for larger redshifts, thus de- 
creasing the detectability of objects located at large redshifts. 
Therefore, most of the observed objects were chosen with re- 
gard to their low redshift (i.e. z < 0.2). 

All observations were carried out in the so called 'wob- 
ble' mode (Aharoni an et al J 19971) . tracking the telescopes with 
an offset of 0.5° in declination with respect to the object po- 
sition, allowing for simultaneous on- and off-source (back- 
ground) observations. For consecutive runs with a duration of 
20 min, the offset sign is reversed in order to avoid system- 
atic effects due to acceptance inhomogeneities in the field of 
vi ew. The background is estimated (similar to the method used 
in lAharonian et alJ (l2003b)) using a ring segment concentric to 
the camera center at the same radial distance to the camera cen- 
ter as the on-source region (i.e. 0.5°). A segment with opening 
angle rf - 70° is excluded from the background region in or- 
der to avoid possible contamination from the on-source region. 
This method makes sure that on- and off-source measurements 
are both taken with identical radial camera acceptances and al- 
lows at the same time large background statistics. 

Two a priori cuts on the system trigger rate are applied to 
each run in order to exclude runs taken under bad weather con- 
ditions and to reduce systematic effects in the determination of 
the excess rate. In a first step a minimum trigger rate of 7 Hz is 
required. This cut excludes data taken under the worst weather 
conditions. In a second step, an expected rate, depending on 
the hardware settings and the zenith angle of the observations 
is calculated from the parameters of a fit to all data of one pe- 
riod with constant hardware settings (all runs with trigger rate 
lower than 7 Hz are excluded from this fit). Runs with rates be- 
low 80 % of the expected rate are rejected. Additionally, runs 
with technical problems are excluded. After the application of 
the above run-selection criteria the total clean data set amounts 
to an observation time of 1017 hours. In Table0all observed 
objects are listed with J2000 coordinates, redshift and object 
type ordered by ascending redshift. Additionally the observa- 
tion time spent on each object as well as the results of the anal- 
ysis described below are listed. 

4. Data Analysis 

The Cherenkov light generated by an air shower initiated by a 
primary y-ray or hadronic particle is seen as an elliptical image 
in each triggered camera. Since each telescope has a different 
viewing angle relative to the shower axis a complete geometri- 
cal reconstruction of the air shower is possible with an image 
analysis of at least two telescopes. 

Before the reconstruction of direction and shower core po- 
sition the following cuts are applied. Reconstructed images 
with more than 15 defective camera pixels are rejected. A min- 
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imum amount of light (size) of 40 photo-electrons (ph.e.) is 
required in an image. Images with a distance of the center of 
gravity to the camera center of more than 78 % of the camera 
radius are rejected in order to avoid truncation by the camera 
border. 



After application of the above image selection criteria at 
least three remaining images are required in this analysis for 
the reconstruction of the direction and the core impact posi- 
tion of an event. This improves the quality of the reconstruc- 
tion, the angular resolution and the separation between y-ray 
and cosmic ray (hadronic background) induced air showers. 
The stereoscopic technique allows for an event-by-event re- 
construction of the direction of the primary particle. Since the 
shape of the elliptical images also depends on the shower core 
position, the reconstruction of the shower impact parameter 
for each telescope provides a means of scaling the individual 
widths of the elliptical shower images of each telescope with 
expected widths for y-ray induced shower images from Monte- 
Carlo simulations. The mean of the scaled widths is called 
mscw-parameter and provides a very good y-hadron separa- 
tion. This is described in detail in Konopelko et al. ( 1999a). 
The optimum cut value for a y-ray signal search is found to 
be mscw =1.1. 



For the reconstruction of the direction of the primary par- 
ticle, algorithm # 3 from ' Hofmann et all ( 1 1999 11) is used. The 
angle A0 = |0o - ©rl between the object direction ©o and the 
reconstructed shower axis 0, is called the angular distance. In 
case of a signal from a source with point-like emission, the dis- 
tribution of the squared angular distance A0^ is expected to 
accumulate entries at small values starting from 0, i.e. the sig- 
nal region. The extension and shape of the signal accumulation 
reflects the angular resolution of the system and depends on the 
telescope multiplicity, the zenith angle and the hardware setup 
of the telescope system. Therefore, the cut on the angular dis- 
tance also depends on the parameters mentioned above. Using 
data of the well known Crab Nebula, the cut on A©- is thus 
optimized individually for different hardware setups, multiplic- 
ities and zenith angle intervals. This method takes the depen- 
dencies of the angular resolution described above into account 
and leads to results consistent with earlier analyses. Typical 
values of the A0^-cut are 0.008 deg^ for events reconstructed 
with five triggered telescopes (having the best angular resolu- 
tion) to 0.015 deg^ for 3-telescope events. Similarly, the cut on 
the core impact position slightly depends on the zenith angle 
(ZA) of the observation. The optimum values found for this cut 
are 200 m (low ZA), 400 m (medium ZA) and 600 m (high ZA). 
In TableQlall selection criteria and cuts are summarized. 



Different cuts on A©^ imply different solid angle ratios of 
on- and off-source region a - Qon/^oS (a-factor) for each sub- 
set. Therefore the significance of an excess is calculated using 
a formula based upon the hkehhood Eq. 17 of .Li & Ma C1983.) 





rate 


> 


7Hz 


run selection 


rate deviation 
technical problems 


< 


20% 




# of defective pixels 


< 


15 


image selection 


image size 


> 


40 ph.e. 




distance 


< 


0.78 




telescope multiplicity 


> 


3 


event selection 


core distance 
mscw 


< 
< 


f(subset) 
1.1 




A0= 


< 


f(subset) 



Table 1. Selection criteria of the analysis chain. The cuts were 
optimized individually for all data subsets using data of the 
well known Crab Nebula (see text). The distance is measured 
from the center of gravity of the image to the camera center 
The entry 'f(subset)' indicates that the cut depends on the data 
subset. 



but generalized for data subsets with different a-factors': 

S = V2 X f y N® ln( ] 

^4^ °" ly„.^(N«+N«)i 



. y N« ln( \ 



"off 



1/2 



The variability of each object is investigated using the 
Kolrn ogorov and the Prahl test (Kolmogorov 1933; Prahl, 
Il999l) . Both tests result in a significance for burst-like be- 
haviour, given a time sequence of events. The Prahl test is espe- 
cially sensitive to burst-like behaviour with a small duty cycle. 

For each object a Crab Nebula y-rate as expected for iden- 
tical observational conditions (zenith angle, hardware setup) is 
calculated from data. These expected rates are used to com- 
pute flux valu es and upper limits on the integral flux following 
lHelenel(ll983h . 

5. Results 

A distribution of the significances for steady state emission 
(DC) of all analyzed objects is shown in Figure ^ The distri- 
bution follows a Gaussian distribution of mean zero and stan- 
dard deviation one, as expected in case of a pure background 
sample, with exceptions from lES 1959-H650, H 1426+428, 
IBS 2344-H514 and M87. In Table H a list of all objects ana- 
lyzed in this work ordered by ascending redshift is shown along 
with their observation time and upper limits on the integral 
flux respectively flux values for the most significant objects. 
The values of the calculated upper limits are found to lie be- 
tween 2 % and 60 % of the Crab Nebula flux. Comparable up- 
per limits for BL Lac objects were found with observations of 
the Whipple telescope, ranging from 6 % to 100% of the Cra b 
Nebula flux faoran et al J 120031: Ide la Calle Perez et allEool . 



' In Li & Ma ( 1983) the significance is derived from the ratio of 
the conditional probabilities for 'background assumption' and 'signal 
assumption'. Substituting both assumptions with a sum over data sub- 
sets with different a-factors and a straight forward calculation leads 
to the above generalized formula. 
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HEGRA data, 1996-2002 
54 AGN 



(1) 1ES 1959+650 

(2) H 1426+428 

(3) 1ES 2344+514 

(4) M 87 



43 2 



Table 2. Number of on- and off-source events and significance 
S for the data sets from the years 1997 (PI), 1998 (P2) and 2002 
(P3) on lES 2344+514. Note that the major part of the excess 
is accumulated in 1998. 



-5 5 10 15 20 25 

DC significance S [a] 

Fig. 1. Distribution of significances for all 54 objects analyzed 
in this work. A Gaussian fit between -2.5 and -1-2.5 cr matches 
the core distribution very well. The objects lES 1959H-650 and 
H 1426+428 show a clear deviation from the background ex- 
pectation. Further two objects, lES 2344+514 and M87 also 
show a deviation from the background expectation on the 4 cr 
level and thus evidence for the emission of TeV y-radiation. 

The statistical tests on burst like variability of all objects yield 
no positive results except for IBS 1959+650 where a high sta- 
tistical significance (> 16 cr) for burst-like behaviour reflects 
the strong obvious flaring activity in May and July 2002. 

The X-ray selected BL Lac o bject IBS 1959+650 
jBlvis et al J h 992: 'Sch achter et al Jll 993h was first reported as 
a TeV y-ray emitter by the Seven Te lescope Array gro up in 
1999 with a DC significance of 3.9 cr jNishivama et all 1 999). 
HBGRA lACT system observations were carried out from 
July to September 2000, May to October 2001 and May to 

September 2002. The HBGR A results on IBS 1959+6 50 were 

I \ II — • — n 

published in detail elsewhere (Aharo nian et al.l '2003c). the re- 
sults presented in this analysis show a level of activity ranging 
from 0.06 to 2.9 Crab units and are consistent with the earlier 
analysis. 

H 1426+428 was re ported to have a synchrot ron peak lying 
near or above 100 keV jCostamante et alJl200ll) thus qualify- 
ing the object as an extreme synchrotron blazar Detections in 
the TeV energy regime fro m this object have been reported by 
the Whippl e collaboration (iHoran et alJ2002l) the HBGRA col- 
laboration ( Aharonian et al. 2002b) and the CAT collaboration 
tD j annati- Atai cl al. 2002 ) . The observations of H 1426+428 
in the years 1999, 2000 and 2002 r esulted in an excess on the 
7.5 cr level ('Aharonian et alJ2003al) . consistent with the analy- 
sis presented here (6.6 cr). 

IBS 2344+514 was one of the first BL Lac type objects 
to be reported as an extreme synchrotron blazar wit h syn- 
chrot ron peak energy reaching up to 100 keV ( Giommi et alJ 



[). The first TeV detect ion of this object wa s reported by 
the Whipple group in 1998 ('Catanese et al."l998'). With an av- 
erage flux of 1 1 % of the Crab Nebula flux in 1998 and a higher 
flux level of 63 % of the Crab flux in one night of observa- 
tions (6 cr), the object has shown clear evidence for a variable 
flux in the Whipple data. The results of the HBGRA obser- 



observation periods 


time 


Non 




S 




[h] 


# 


# 


[0-] 


Pi Oct - Dec 1997 


15.0 


54 


52 


0.3 


P2 Aug - Nov 1998 


41.8 


128 


84 


4.3 


P3 Sep 2002 


15.7 


53 


35 


2.6 


2 


72.5 


235 


171 


4.4 



vations on IBS 2344+514 of the year 19 97 and 1998 were 
first reported by iKonopelko et al.l d 1999 b') with a DC signif- 
icance of 3.3 cr. Further observations have been carried out 
since the above publication. The analysis presented here in- 
cludes the complete dataset and results in an excess of 64 + 15 
photons (Non - 235, <Noff> - 171) with a significance of 
4.4 cr. The data set of IBS 2344+514 can be spHt into three 
independent observation periods. The first period PI ranges 
from October to December 1997, the second period P2 from 
August to November 1998. PI and P2 are separated by a pe- 
riod of non observability of the object from the HBGRA site. 
Additional observations have been carried out in September 
2002 (P3). The data subset PI shows no evidence for a TeV 
y-ray signal, with a DC significance of 0.3 cr whereas the sec- 
ond observation period P2 yields a significance of 4.3 cr. In 
the last observation period P3 an excess on the 2.6 cr level 
is found. Tests for burst-like behaviour do not yield statisti- 
cally significant results. In Table|2lthe number of on- and off- 
source events as well as the corresponding significances are 
listed for the different data subsamples. In Figure|2the distri- 
butions of the reconstructed directions for the complete data set 
(P1+P2+P3) and the data set with the highest significance (P2) 
of IBS 2344+514 are shown. The observed excess results in a 
flux of <D(B>0.97 TeV) = (0.60 + 0. 19) ■ lO^'^ photons cm^^s-i , 
corresponding to (3.3 + 1.0)% of the Crab Nebula flux. The 
errors on the flux level are dominated by the statistics of the 
measurements. 

In contrast to the abovementioned 4 objects the jet of the 
giant radio galaxy M87 is not aligned to our line of sight which 
makes it the only non BL Lac type object among the 4 most sig- 
nificant objects of this data sample. The observed excess from 
the radio galaxy M87 results in a significance of 3.9 cr in the 
present analysis and is consistent to the results of a detailed 
analysis of the M87 data, which yield 4.1 cr a nd were presented 
in a dedicated paper dAharonian et alJ20 03b'). After further im- 
provements this analysis now yields 4.7 cr ( Getting et al. 

Several objects that were reported to be sources of TeV y- 
rays by other groups are included in this data set. Among these, 
the objects 3C066A (Neshporet al. 1998) and PK S 2155-304 
dChadwick et al.l ll 999^ IPi annati-Atai et all l2003h were only 
observed for a very short time. No excess was found in the 
HBGRA data of these two objects. A weak excess on the 3 cr 
level is found in BL-Lacertae which was observed for 29 hours, 
resulting in a 99 % C.L. upper limit of the order of 28% of the 
Crab Nebula flux (see Table^}. This object was reported to be a 
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1ES 2344+514 



HEGRA data, 1997-2002 




0.01 0.02 0.03 0.04 0.05 0.06 



squared angular distance Ae^ [deg^] 

Fig. 2. Distributions of reconstructed squared angular distances 
to the object direction (see text) of the lES 2344+514 data. 
The distribution of the on-source events is represented by the 
data points. The scaled background (off) is shown as a shaded 
histogram. The upper figure represents the total data set taken 
between 1997 and 2002 whereas the lower figure shows the 
1998 data alone. As can be seen from this figure the excess 
accumulates essentially in the year 1998. 

TeV y -ray emitter by the Crimean Observatory jNeshpor et alJ 
6. Discussion 

Predictions for y-ray fluxes in the GeV/TeV energy reg ime 
have been made for several AG N bv lStecker et alJ ll 19961) and 
ICostamante & Ghisellini' ('2002'). 

In Table|3the observed integral fluxes and flux upper limits 
of several objects of the present data sarn ple are compared to 
predictions made by Stecke r et al ] 11996^, based upon simple 
scaling arguments, taking into account the Einstein Slew survey 
sample of BL Lacs. The authors argue that only high frequency 
peaked BL Lac objects are potential sources of extragalactic 
TeV y-radiation. For all objects listed in Table |3l the derived 
HEGRA upper limits on the integral flux were calculated for 
an energy threshold of 1 TeV and for different spectral indices, 
assuming a power law energy spectrum. The observed fluxes of 



Table 3. Comparison of upper limits and fluxes (O) derived 
in this work with predictions (o„k,jci) made by Stecker et alJ 
(1996). All fluxes are given in units of 10"'^ photons cm^^s"'. 
The 99 % C.L. upper limits were extrapolated to a fixed energy 
threshold of 1 TeV assuming a power law energy spectrum for 
3 differential spectral indices a = 2.0, 2.5, 3.0. 



object z 


Dmodcl/ 10-'- 




o/io-'- 






ycm-- s"' 




y cm"- 


1 






spectral index a 






2.0 


2.5 


3.0 


IZwl87 0.055 


0.59 


<1.56 


<1.51 


<1.47 


IBS 2321+419 0.059 


0.14 


<0.60 


<0.56 


<0.53 


lES 1741 + 196 0.083 


0.35 


<1.33 


<1.29 


<1.25 


PKS 2155-304 0.116 


0.88 


<1.60 


<3.83 


<9.16 


lES 1118+424 0.124 


0.18 


<4.18 


<4.12 


<4.06 


1ES0145+138 0.125 


0.26 


<1.19 


<1.11 


<1.04 


lES 1212+078 0.130 


0.03 


<2.98 


<2.86 


<2.74 


lES 0229+200 0.139 


0.11 


<2.99 


<2.87 


<2.75 


lES 1255+244 0.140 


0.34 


<2.03 


<1.97 


<1.91 


lES 0323+022 0.147 


0.15 


<0.71 


<0.71 


<0.71 


lES 1440+122 0.162 


0.03 


<1.37 


<1.31 


<1.26 


lES 0347-121 0.185 


0.08 


<7.50 


<9.07 


<10.96 


lES 0927+500 0.186 


0.02 


<1.02 


<0.98 


<0.95 


lES 2344+5 14 0.044 


0.80 


0.58 


0.57 


0.56 


lES 1959+650 0.047 


2.30 


1.0-51 
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lES 1959+650 in its low state and lES 2344+514 are close to 
the predicted values. The calculated HEGRA upper limits for 
all other objects in this list exceed the predictions. For some 
objects (e.g. lES 0927+500, lES 1440+122) the predicted flux 
levels are a factor of 50 to 100 lower than the upper Umits from 
HEGRA. The level of sensitivity necessary to detect such low 
fluxes within a reasonable time is beyond the capabilities of the 
HEGRA Cherenkov telescopes. Such low fluxes could only be 
detected by the experiments of the next generation which have 
a higher sensitivity and a lower energy threshold. 

In Table |3 the HEGRA integral flux upper limits de- 
rived in this w ork are compared to predictions made by 
ICostamante & Ghisellini (2002). Two different model predic- 
tions are given. The first number is taken from a parametriza- 
tion of the spectral energy distribution (SED) origin ally intro- 
duced by Fossati et al.. (1998.) and modified by .Donato et alJ 
(1200 ll) and ICostamante & Ghiseliiiiil (l2002l) (hereafter FDC). 
The second number is calculated using an SSC model from 
ICostamante & Ghisellini C200Z) (hereafte r CG). Compared to 
the earlier work of IStecker et al I (Il996l) . CG introduce the 
new requirement of strong radio emission for a TeV candidate 
source (arguing that a strong radio emission is a good indicator 
for non-thermal low energy emission producing seed photons). 
Furthermore, several other BL Lac samples in addition to the 
Einstein Slew survey sample were taken into account. The FDC 
parametrization is rather suitable for predictions of high state 
TeV fluxes while the SSC model predictions, designed to fit 
the known synchrotron part of the SED, are more appropriate 
for a quiescent state of the TeV source candidate. Additionally, 
one has to note that along with other uncertainities the absorp- 
tion of TeV photons by the extragalactic background radiation 
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Tab le 4. Comparison of flux predictions by 
ICos tamante & Ghisellini (2002) (*m„dd) with the results 
derived in this work. All fluxes are given in units of 
10 '^photons cm^^s '. For the extrapolation of the upper 
limits given in Table Q a power law for 3 differential spectral 
indices (a - 2.0, 2.5, 3.0) were assumed and the energy 
threshold was fixed at Ethr = 1 TeV. For lES 0647+250 a red- 
shift of 0.200 was a ssumed Two numbers for predicti ons above 
ITeV are given bv ICostamante & Ghisell ini (2002): the first 
num ber is obtained from a parametrization of the SED adapted 
fro mlFossati et alJ( 1998) and modified by Donato et al. ( 2001) 
and ICostamante et alJ ( 120011) (FDC). The second number 
resu lts from a homogene ous one -zone SSC model described in 
ICost amante & Ghisellinil (Eool) (CG). 



object 


z 








<s>/n)-'- 








7 cm 


-2 S-' 


7 cm - s 








FDC 


CG 


spectral index a 










2.0 


2.5 


3.0 


H 1722+119 


0.018 


35.2 


0.01 


<3.84 


<3.62 < 


3.41 


lES 1741+196 0.083 


8.4 


0.1 


<1.33 


<1.29< 


1.25 


lES 0806+524 


0.138 


2.7 




<4.63 


<4.84 < 


5.05 


lES 0229+200 


0.139 


2.1 


0.04 


<2.99 


<2.87 < 


2.75 


RBS 0958 


0.139 


2.8 




<5.30 


<4.88 < 


4.50 


lES 0323+022 0.147 


1.8 




<0.71 


<0.71 < 


0.71 


lES 1440+122 0.162 


2.0 


0.1 


<1.37 


<1.31 < 


1.26 


PG 1218+304 


0.182 


1.5 




<2.24 


<2.06 < 


1.88 


lES 0647+250 


0.200 


1.2 




<2.61 


<2.31 < 


2.04 


lES 1011+496 


0.200 


0.2 




<1.84 


<1.85 < 


1.87 


lES 0120+340 


0.272 


0.6 




<0.72 


<0.66 < 


0.60 


lES 1959+650 


0.047 


17.4 - 


1.0-51.9 



field was not accounted for in these models (CG). In addition 
to the flux predictions given in CG, we have used the FDC 
parametrization to calculate flux predictions for several objects 
not included in the list of CG. We have additionally included 
the effect of the absorption of TeV photons by pair production 
with the extragalactic background light (EBL). For this purpose 
we have used a model parametrization of the spectral energy 
distribution of the EB L (also us e d in A haronian et al. (2002b), 
model 1) adopted from lPrimackl J2OO1]) and designed to be con- 
sistent with our observations of known TeV Blazars. The lu- 
minosity distance was calculated following Ue^Li Pen ( 1996). 
These predictions are compared to upper limits calculated in 
this work in Table|5] In the case of the FDC parametrization, the 
predicted flux levels exceed the derived upper limits for those 
objects printed in boldface in Table|3 and Table|5l Assuming 
these predictions to hold true for a high state of activity, it can 
be concluded that these objects were not in a flaring state during 
the HEGRA observations. However, if we take into account the 
absorption by the EBL (column 4, Table|5} this only remains 
valid for four objects. 

For lES 1959+650 the observed flux level during the high- 
est state of emission in the HEGRA data is found to ex- 
ceed the predicted value by a factor of 3 (only FDC) to 12 
(FDC+EBL). Given the high variability of the object at this 
time, prediction and observation can easily be accomodated. 
The observed flux from H 1426+428 is roughly a factor of 3 



Table 5. Comparison of integral flux predictions above 1 TeV 
using an own implementation of the FDC parametrization 
and including the absorption by the extragalactic background 
light (Omodd) with the 99 % C.L. flux upper limits (O) derived 
in this work. For the objects lES 1959+650, H 1426+428 and 
lES 2344+514 fluxes are given. Objects printed in boldface 
show upper limits below the predicted values. In case of 
lES 0647+250 a redshift of 0.200 was assumed in order to be 
able to calculate the EBL absorption. 



object 


z 


•t'modcl 


/lo-i- 




o/io-'- 








7 cm 


-2 S-' 


7 


cm-^ 5- 


1 






FDC 


FDC 


spectral index a 








+EBL 


2.0 


2.5 


3.0 


H 1722+119 


0.018 


36.98 


20.53 


<3.84 


<3.62 


<3.41 


lES 2344+514 


0.044 


18.58 


4.72 


= 0.58 


= 0.57 


= 0.56 


Mkn 180 


0.046 


19.00 


4.56 


<1.64 


<2.00 


<2.45 


lES 1959+650 


0.047 


18.01 


4.20 


= ( 


1.0-51 


.9) 


3C371 


0.050 


12.41 


2.80 


<2.51 


<3.09 


<3.81 


I Zw 187 


0.055 


12.40 


2.25 


<1.56 


<1.51 


<1.47 


lES 2321+419 


0.059 


4.85 


0.75 


<0.60 


<0.56 


<0.53 


BL-Lacertae 


0.069 


1.66 


0.22 


<4.51 


<4.73 


<4.96 


lES 1741 + 196 


0.083 


7.10 


0.55 


<1.33 


<1.29 


<1.25 


PKS 2155-304 


0.116 


2.40 


0.07 


<1.60 


<3.83 


<9.16 


lES 1118+424 


0.124 


2.36 


0.05 


<4.18 


<4.12 


<4.06 


lES 0145+13.8 


0.125 


1.16 


0.02 


<1.19 


<1.11 


<1.04 


IH 0658+595 


0.125 


2.49 


0.05 


<0.98 


<1.02 


<1.06 


H 1426+428 


0.129 


2.32 


0.04 




= 0.69 




lES 1212+078 


0.130 


2.74 


0.05 


<2.98 


<2.86 


<2.74 


lES 0806+524 


0.138 


1.85 


0.03 


<4.63 


<4.84 


<5.05 


lES 0229+200 


0.139 


2.25 


0.03 


<2.99 


<2.87 


<2.75 


RBS 0958 


0.139 


2.42 


0.03 


<5.30 


<4.88 


<4.50 


lES 1255+244 


0.140 


1.12 


0.01 


<2.03 


<1.97 


<1.91 


lES 0323+022 


0.147 


1.97 


0.02 


<0.71 


<0.71 


<0.71 


OQ530 


0.152 


0.23 


3 10-^^ 


<1.58 


<1.67 


<1.77 


lES 1440+122 


0.162 


1.69 


0.01 


<1.37 


<1.31 


<1.26 


PKS 0829+046 


0.180 


0.03 


210^^ 


<0.96 


<0.96 


<0.96 


PG 1218+304 


0.182 


1.24 


5 10-^ 


<2.24 


<2.06 


<1.88 


lES 0347-121 


0.185 


0.83 


3 10-^ 


<7.50 


<9.07 


<10.96 


lES 0927+500 


0.186 


1.06 


3 10-^ 


<1.02 


<0.98 


<0.95 


PKS 2254+074 


0.190 


0.06 


2 10-* 


<0.89 


<0.85 


<0.80 


MS 0317+1834 


0.190 


1.00 


3 10-^ 


<2.37 


<2.12 


<1.89 


lES 0647+250 


0.200 


0.83 


210^^ 


<2.61 


<2.31 


<2.04 


lES 1011+496 


0.200 


0.29 


7 10-" 


<1.84 


<1.85 


<1.87 


lES 0120+340 


0.272 


0.37 


810-5 


<0.72 


<0.66 


<0.60 


2E 0414+0057 


0.287 


0.22 


3 10-5 


<2.18 


<2.19 


<2.20 


S50716+714 


0.300 


0.21 


210-5 


<4.95 


<6.22 


<7.81 


3C 66A 


0.444 


0.11 


1 10-' 


<3.29 


<3.03 


<2.80 


PKS 0219-164 


0.698 


0.04 


5 10-'2 


<3.29 


<4.39 


<5.86 



below the value predicted by the FDC parametrisation alone. 
But including the absorption by the EBL the predicted value 
is much lower than the observed flux (FDC+EBL). The ob- 
served flux level from lES 2344+514 is lower than the pre- 
dicted value by a factor of 30 (only FDC) resp. 8 (FDC+EBL). 
However, this object has shown flux levels in the earlier 
Whipple data which exceeded the flux observed by HEGRA 
by a factor of 20, which shows that during HEGRA observa- 
tions lES 2344+514 was indeed not in a flaring state. The SSC 
model of ICostamante & GhiseUinL (,200Z) predicts flux values 
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well below the observed upper limits. Most of these predicted 
values would only be detectable with much longer exposure 
times or with the next generation of Cherenkov telescopes. 

7. Summary 

A data set of 54 Active Galactic Nuclei observed with the 
HEGRA lACT system between 1996 and 2002 has been anal- 
ysed. The two objects lES 1959+650 and H 1426+428 are 
meanwhile well established sources of TeV y-radiation. The 
evidence for TeV y-rays from IBS 2344+514 is a confirmation 
of the detection of this object by the Whipple collaboration. 
The detection of TeV y-rays from M87 would be, if confirmed, 
the first detection of photons in the TeV energy regime from an 
AGN of an object not commonly classified as a BL Lac object. 
Upper limits have been derived for all other 50 objects. Table|6l 
summarizes the results for the most significant excesses seen in 
this data set. Comparisons with diff'erent model predictions in- 
dicate that a higher sensitivity is needed to be able to constrain 
the models. 

A further step towards understanding the involved accel- 
eration mechanisms and of the AGN class as a whole as well 
as of the absorption by the extragalactic background radiation 
field is expected from further observations of AGN over a wide 
range of redshifts and especially BL Lac type objects with the 
next generation of Cherenkov telescopes (partly already in op- 
eration) and with future instruments. With the reduction of the 
energy threshold towards lOOGeV the eff'ect of the absorption 
by the extragalactic background light will become less impor- 
tant for objects located at low redshift. Thus, the uncertainties 
in the interpretation of observations induced by this effect will 
be less important. 
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Table 6. This Table summarizes the resuhs of the present analysis for the most significant objects from the sample of 54 AGN 
observed with the system of stereoscopic Cherenkov telescopes. The number of on- and off-events as well as significance S, 
fluxes in units of the Crab Nebula flux (F) and in units of 10 '^photons cm^^s ' (O) above the energy threshold Ethr of the 
observation are given. 
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[10"'^ photons cm^^s"'] 
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0.6 - 30.5 
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0.91 


0.03 


0.8 
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Table 7. List of all objects of the HEGRA AGN data sampl e. The J2000 co ordinates are given as well as redshift and object type 
(following Stickel et al. ( 1994); Padovani & Giommil d 1 995l) :IVeron-Cettv & Veronl JioOl)). The results of the analysis presented 
here are summarized. In case of the objects deviating significantly from the background expectation the flux is given in units of 
the Crab Nebula flux. For all other objects upper limits on the integral flux in units of the Crab Nebula flux (F^j?*) and in units of 



10 '^photons cm ^s 



(O^'fO are given. 



object 


Or 


6 


z 


type 


T 


Ethr 


S 




o?,7/io- 




(hh mm ss) 


(dd mm ss) 






[h] 


[TeV] 


[tr] 


[Crah] 


[y cm~^ s 


lES 0647+250 


06 50 46.6 


+25 03 00 


— 


BL 


4.1 


0.78 


0.5 


0.13 


3.35 


MG 0509+0541 


05 09 26.0 


+05 41 35 


— 


BL 


15.8 


0.96 


1.3 


0.11 


1.92 


M87 


12 30 49.4 


+ 12 23 28 


0.004 


Fl 


70.0 


0.88 


3.9 






NGC315 


00 57 48.9 


+30 21 08 


0.016 


Fl/2 


14.6 


0.86 


-0.3 


0.05 


1.03 


NGC 1275 


03 19 48.2 


+41 30 42 


0.018 


Fl 


87.6 


0.85 


-0.3 


0.03 


0.68 


H 1722+119 


17 25 04.5 


+ 11 52 15 


0.018 


BL 


5.1 


0.89 


1.7 


0.21 


4.31 


PKS 2201+04 


22 04 17.7 


+04 40 03 


0.028 


SI 


17.8 


0.95 


1.7 


0.08 


1.40 


VZw331 


03 13 57.0 


+41 15 37 


0.029 


BL 


4.1 


0.87 


-0.2 


0.09 


1.93 


NGC 1054 


02 42 15.0 


+ 18 13 00 


0.032 


G 


57.9 


0.86 


-1.7 


0.02 


0.37 


3C 120 


04 33 12.0 


+05 21 15 


0.033 


Fl 


25.4 


0.93 


-0.7 


0.05 


0.86 


NGC 4151 


12 10 32.7 


+39 24 19 


0.033 


S1.5 


7.0 


0.79 


-0.4 


0.07 


1.79 


UGC01651 


02 09 38.5 


+35 47 51 


0.037 


G 


14.3 


0.79 


1.3 


0.07 


1.62 


UGC03927 


07 37 30.0 


+59 41 03 


0.041 


F2 


6.3 


1.09 


-2.4 


0.09 


1.32 


lES 2344+514 


23 47 04.9 


+51 42 17 


0.044 


BL 


72.5 


0.97 


4.4 






MknOlSO 


11 36 26.4 


+70 09 27 


0.046 


BL 


9.8 


1.50 


-0.6 


0.12 


1.09 


lES 1959+650 


19 59 59.9 


+65 08 54 


0.047 


BL 


163.7 


1.32 


28.0 






3C 371.0 


18 06 50.7 


+69 49 28 


0.050 


BL 


5.4 


1.52 


-0.4 


0.19 


1.65 


4C +37.11 


04 05 49.3 


+38 03 32 


0.054 


S 


6.7 


0.80 


-2.0 


0.05 


1.17 


IZw 187 


17 28 18.6 


+50 13 10 


0.055 


BL 


16.0 


0.94 


1.9 


0.09 


1.66 


Cyg-A (3C 405.0) 


19 59 28.5 


+40 44 02 


0.057 


F2 


59.0 


0.91 


-0.2 


0.03 


0.64 


lES 2321+419 


23 23 52.5 


+42 10 55 


0.059 


BL 


22.3 


0.89 


-1.6 


0.03 


0.67 


3C 192.0 


08 05 35.0 


+24 09 50 


0.060 


F2 


2.9 


0.93 


0.3 


0.20 


3.78 


4C+31.04 


01 19 35.0 


+32 10 50 


0.060 


FR 


3.0 


0.76 


-0.3 


0.14 


3.83 


BL-Lacertae 


22 02 43.3 


+42 16 40 


0.069 


BL 


26.7 


1.10 


3.0 


0.28 


4.10 


lES 1741 + 196 


17 43 57.8 


+ 19 35 09 


0.083 


BL 


10.2 


0.94 


0.3 


0.07 


1.41 


4C+01.13 


05 13 52.5 


+01 57 10 


0.084 


F2 


7.7 


1.01 


-0.2 


0.10 


1.73 


PKS 2155-304 


21 58 52.0 


-30 13 32 


0.116 


BL 


1.8 


5.72 


0.0 


0.27 


0.28 


lES 1118+424 


11 20 48.1 


+42 12 12 


0.124 


BL 


2.0 


0.97 


0.3 


0.24 


4.31 


lES 0145+13.8 


01 48 29.8 


+ 14 02 19 


0.125 


BL 


3.2 


0.87 


1.1 


0.06 


1.37 


IH 0658+595 


07 10 30.1 


+59 08 20 


0.125 


BL 


33.7 


1.08 


-0.4 


0.06 


0.91 


H 1426+428 


14 28 32.5 


+42 40 25 


0.129 


BL 


258.5 


0.91 


6.6 






3C197.1 


08 21 32.6 


+47 02 46 


0.130 


QSO 


15.0 


0.96 


-0.4 


0.05 


0.86 


lES 1212+078 


12 15 11.2 


+07 32 02 


0.130 


BL 


2.4 


0.92 


-0.6 


0.17 


3.24 


lES 0806+524 


08 09 49.2 


+52 18 58 


0.138 


BL 


1.0 


1.09 


-0.1 


0.29 


4.25 


lES 0229+200 


02 32 48.7 


+20 17 17 


0.139 


BL 


3.0 


0.92 


1.0 


0.17 


3.25 


RBS 0958 


11 17 06.3 


+20 14 06 


0.139 


BL 


3.8 


0.85 


2.7 


0.28 


6.23 


lES 1255+244 


12 57 32.0 


+24 12 39 


0.140 


BL 


5.9 


0.94 


0.1 


0.12 


2.16 


MS1019.0+5139 


10 22 11.0 


+51 24 00 


0.141 


S 


17.5 


0.92 


0.1 


0.07 


1.35 


lES 0323+022 


03 26 13.9 


+02 25 14 


0.147 


BL 


14.3 


1.00 


-1.5 


0.04 


0.71 


GQ530 


14 19 46.6 


+54 23 14 


0.152 


BL 


9.4 


1.12 


0.4 


0.10 


1.41 


3C 273.0 


12 29 06.7 


+02 03 08 


0.158 


FR 


12.2 


1.15 


-0.3 


0.09 


1.25 


lES 1440+122 


14 42 48.4 


+ 12 00 39 


0.162 


BL 


13.1 


0.92 


-0.9 


0.08 


1.49 


PKS 0829+046 


08 31 48.9 


+04 29 39 


0.180 


BL 


18.0 


1.00 


0.5 


0.06 


0.96 


PG 1218+304 


12 21 22.0 


+30 10 37 


0.182 


BL 


3.9 


0.84 


-0.3 


0.12 


2.67 


lES 0347-121 


03 49 23.0 


-11 59 26 


0.185 


BL 


1.9 


1.46 


1.8 


0.56 


5.14 


lES 0927+500 


09 30 37.6 


+49 50 24 


0.186 


BL 


13.3 


0.94 


0.2 


0.06 


1.08 


PKS 2254+074 


22 57 17.3 


+07 43 12 


0.190 


BL 


16.3 


0.90 


-0.5 


0.05 


0.99 


MS0317.0+1834 


03 19 51.9 


+ 18 45 35 


0.190 


BL 


2.7 


0.80 


-0.5 


0.12 


2.96 


lES 1011+496 


10 15 04.2 


+49 26 00 


0.200 


BL 


2.0 


1.02 


-1.3 


0.11 


1.80 


lES 0120+340 


01 23 08.9 


+34 20 50 


0.272 


BL 


18.9 


0.83 


-1.2 


0.04 


0.87 


2E 0414+0057 


04 16 52.5 


+01 05 23 


0.287 


BL 


4.5 


1.01 


0.6 


0.13 


2.16 


S5 0716+714 


07 21 53.4 


+71 20 36 


0.300 


BL 


1.7 


1.58 


0.7 


0.38 


3.13 


3C 066A 


02 22 39.6 


+43 02 07 


0.444 


BL 


1.3 


0.85 


-0.2 


0.17 


3.87 


PKS 0219-164 


02 22 01.0 


-16 15 16 


0.698 


BL 


1.7 


1.78 


-1.7 


0.27 


1.85 



F 

[Crab] 



0.04 



0.03 
0.06-2.9 



0.03 



